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Abstract. Radio telescopes with ofF-axis feeds, such as the (E)VLA, suffer from "beam squint" in which the 
two orthogonal circular polarizations sampled have different pointing centers on the sky. Its effects are weak 
near the beam center but become increasingly important towards the edge of the antenna power pattern where 
gains in the two polarizations at a given sky position are significantly different. This effect has limited VLA 
measurements of circular polarization (Stokes V) and introduced dynamic range limiting, wide-field artifacts in 
images made in Stokes I. We present an adaptation of the visibility-based deconvolution CLEAN method that can 
correct this defect "on the fly" while imaging, correcting as well the associated self-calibration. We present two 
examples of this technique using the procedure "Squint" within the Obit package which allows wide-field imaging 
in Stokes V and reduced artifacts in Stokes I. We discuss the residual errors in these examples as well as a scheme 
for future correction of some of these errors. This technique can be generalized to implement temporally- and 
spatially-variable corrections, such as pointing and cross-polarization leakage errors. 

Key words. Techniques: image processing. Techniques: interferometric 



1. Introduction 

The enhanced sensitivity of the EVLA will allow it to 
study sources much weaker than are currently observed. 
For frequencies at the lower end of the EVLA range this 
will generally require imaging the full primary beam of 
the antenna pattern, or at least those portions contain- 
ing the stronger sources in the field in order to remove 
their sidelobes from the region of interest. In this regime, 
numerous effects of minor importance to observations of 
strong sources near the antenna pointing position must 
be understood and removed in order to allow the EVLA 
to reach its potential. "Beam squint" resulting from the 
off-axis arrangement of the feeds is one of these effects. 

Many galactic objects have significant circular polar- 
ization in their emission. Without correction for the beam 
squint, the (E)VLA has strong off-axis instrumental cir- 
cular polarization (Figure [1]) which masks any intrinsic 
circular polarization. Correction for the beam squint will 
result in a significant increase in the sensitivity of mea- 
surements of circular polarization. 

(E)VLA observations are performed using dual orthog- 
onal circular polarization (henceforth RCP and LCP). It 
has been known for quite some time that the antenna gain 
patterns formed on the sky for these two polarizations are 
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not concentric, having a relative offset of 0.06±0.005 of the 
antenna power FWHM; i n reasonable agreement w ith a 
calculated value of 0.053 (JNapier fc Gustincidll977l ). The 
NVSS survey included careful measurements of the pri- 
mary beam response of the VLA antennas at a frequency 
of 1.4 GHz. The beam squint was found to be 1.71' ±0.02' 
(jCondon et al.l Il998f) corresponding to 0.055 x FWHM 
(Figure [T|). Thus, for perfectly pointed observations, the 
gain in RCP and LCP towards a source located away 
from the pointing position will be modulated by up to 
~ 7% due to the Earth's rotation in approximately oppo- 
site sense below or above the nominal gain depending on 
the location of the source. The effect will be approximately 
equal in magnitude for the RCP and the LCP but of op- 
posite sign and thus will largely cancel when estimating 
Stokes I; as long as both polarizations have been observed 
and preserved in the editing process. However, amplitude 
self-calibration might magnify the effect if corrections are 
determined separately for both polarizations as they will 
tend to counter the modulation of the strongest source 
affected by the beam squint. This can be prevented by 
preaveraging the RCP and LCP data prior to determin- 
ing common amplitude correction factors. 

In this paper, we discuss the beam squint in some 
detail, present an algorithm that corrects its effects 
and apply it to observations made with the VLA. 
We discuss as well the possible extension of this 
technique to a wider class of time and/or direction- 
dependent corrections. All data manipulation discussed 
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Fig. 1. The VLA primary an tenna pattern as m ea- 



1998h 



sured during the NVSS survey (jCondon et al.l 119981) at 
1.4 GHz. The power pattern is shown in gray-scale with 
a scale bar at the top and contours are plotted at the 
±0.02, ±0.04, • • ■ ± 0.16 levels of fractional Stokes V (from 
the center outwards). Negative contours are dashed. 



in this report use the Obit package (jCottonl 120081 ) 
( http://www.cv.nrao.edu/~bcotton/Obit.html I . 



2. The (E)VLA beam squint 

Extensive studies of the polarization properties of "off- 
set parabolas" have been performed since the early 1970s 
with detailed calculations describing the squint and cross- 
polarization propert ies of prime-focus offset parabolas 
(jChu fc TurrinI 119731 ). Subsequent, more refined studies 
determined the polarization properties of a single-reflector 
parabola illuminated by a circularly-polarized feed located 
at the prime focus and whose axis is tilted by an angle 
00 with respect to the axis of the parabola. The result- 
ing antenna primary beam is directed off-axis orthogo- 
nally to the plane that contains the axis of the parabola 
as well as that of the feed. T he beam squint is given by 
( Duan fc Rahmat-Samiilll99ll ): 



sin9oX 
AttF 



(1) 



where 9q is the angle between the feed axis and the axis 
of the parabola, A is the wavelength of the observations 
and F is the focal length of the parabola, (i^ = 9.0 m at 
the (E)VLA) and =F Os is the resulting beam squint of the 
(LCP, RCP) beam. 

The dual reflector offset Cassegrain antenna was solved 
by defining an equivalent prime-focus parabola that is 



equivalent in performance to the wanted dual reflector 
Cassegrain antenna as it has the same co-polarized and 



cross-polarized radiation patterns ( Rusch et al.lll990l ). It 



can be shown that the beam squint of such an antenna 
configuration is given by 



sin{9f} — a)X 



AttF, 



oq 



(2) 



where a is the angle between the axis of the equivalent 
parabola and that of the feed and Op is the angle between 
the feed and the vertex of the off-axis, shaped subreflector 
and Feq is the equivalent focal length of the dual reflector. 
At the (E)VLA antennas, the corres ponding valu es are 



a = 59.8°, Op = 1° and Feq = 59.1 m (|NapieiJll994l) . This 



reference discusses the derivation in some detail and shows 
that the beam squint is thus described by an inverse-sine 
function. However, the small-angle approximation is ap- 
propriate and thus the squint angle at the (E)VLA an- 
tennas scales quasi-linearly with wavelength. This is not 
necessarily the case with the FWHM of the primary beam 
response as this depends on the specific illumination due 
to the feed horns at each band. The description of the 
primary beam response as corresponding to uniform illu- 
mination of a 24.5 m circular aperture (appropriate at the 
VLA at 1.4 GHz given the 25 m antenna diameter and 
the taper applied at the edge of the primary) leads to the 
5.5% ratio discussed above. We have found that this ra- 
tio holds as well for the new (E)VLA 5 GHz feed system 
(Uson & Cotton, in preparation) but do not know at this 
time what the corresponding ratio at other frequencies of 
the (E)VLA will be. We intend to perform the relevant 
measurements in due course. 

3. Beam squint correction technique 

The beam squint can be characterized as a coupled offset 
of the pointing of the two orthogonally polarized beams. 
The antenna will be nominally pointed using the midpoint 
of these two beams giving an (opposite) effective point- 
ing error in each polarization. In the following, the beam 
squint is treated as a polarization-dependent pointing er- 
ror. 

The beams are offset orthogonally to the line connect- 
ing the feed and the antenna center. Due to the alt-az 
mounts of the (E)VLA antennas, this offset will rotate on 
the sky with parallactic angle. Notice that by construction 
of the station pads, all (E)VLA antennas have the same 
parallactic angle at any given time even in the extended 
configurations. 

The effect of beam squint on a given source depends 
on its location in the field of view so, in general, no op- 
eration on the uv data can remove the effects of beam 
squint over the whole field, although an approximation can 
be attempted t hrough construction o f a suitable unitary 
operator (e.g. (JBhatnagar et al.l 120081 )1. However, such a 
procedure encounters problems on fields with significant 
emission from sources located near the nulls of the pri- 
mary beam. The effects of beam squint can be included 
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whenever the instrumental response to a model (e.g., a 
set of CLEAN components) is computed. In the algorithm 
described here, this is done using a discrete Fourier trans- 
form where the instrumental response for each visibility 
is computed for each component and then summed over 
components. The flux density of each component is modi- 
fied by the ratio of the antenna gain in the nominal point- 
ing direction to the actual, beam squint dependent, gain. 
This is done independently for each polarization. As the 
beam squint rotates with parallactic angle, this flux den- 
sity correction is time as well as direction dependent. 

There are two instances where this interferometer 
response model is used. The first is in a vi sibility- 
based f "Schwab-Cotton" . ISchwabI (|l983h . lCottonl |l9991 1 
CLEAN in which a set of compon ents is i nitial ly located 
by a major cycle of the "Clark" , IClarkI (|l980l ) CLEAN 
followed by an accurate model calculation allowing the ef- 
fects of those components to be removed from the visibility 
data and a new residual image derived from the residual 
visibility data. Thus, errors in the initial image due to in- 
adequacies in the initial response model are corrected as 
deconvolution progresses. 

The other use of the instrumental response model is 
in the self-calibration step in which the visibilities are di- 
vided by the corresponding instrumental response to the 
sky model to produce a dataset equivalent to observations 
of a point source with a perfect instrument. Using an ac- 
curate beam-squint model allows the removal of the effects 
of beam squint from the resultant gain solution. 

If an image in Stokes V is desired, then Stokes I must 
be imaged first and the Stokes I model subtracted from the 
data with beam squint corrections applied. This should 
remove the instrumental circular polarization from the 
data. Subsequently, the Stokes V image can be derived 
and deconvolved without further beam squint corrections. 
Instrumental leakage terms (the so-called D-terms) have 
been ignored in the procedure discussed in this paper, but 
could be incorporated if necessary as long as such leakage 
terms were determined with sufficient precision. 



3.1. CLEAN model for beam squint corrections 

In the normal model calculation, a list of CLEAN com- 
ponents is kept and, for each visibility measurement, the 
instrumental response to each component is calculated 
with the total instrumental response obtained by summing 
over all components. To apply the beam squint correction, 
the time- and polarization-dependent antenna gain correc- 
tions are used to adjust the component flux density for the 
effects of beam squint on the antenna power pattern. To 
do this, the antenna voltage gain relative to the nominal 
pointing is determined for the RCP and LCP beams in 
the direction of each CLEAN component and the intcrfcr- 
ometric (power) gain correction is obtained as the product 
of the antenna voltage gain corrections. These gain correc- 
tions are updated whenever the observing parallactic angle 



changes by one degree. Of course, this interval could be 
shortened if necessary. 

A mixed array such as the VLA-I-EVLA adds another 
complication as t he feed placement is diff erent for the two 
types of antennas (ICotton fc Usonll2007ah . In this case, the 
right- and left-hand voltage gain ratios are computed sep- 
arately for VLA and EVLA antennas and the appropriate 
voltage gains are applied to the component flux density 
for each baseline. 

The implementation described above assumes that all 
antennas observe the sky at the same parallactic angle. 
This is indeed the case at the (E)VLA as care was taken 
during the construction of the VLA to build the antenna 
mounts to a common horizon. For arrays such as the 
VLBA, parallactic angle will differ from antenna to an- 
tenna and an implementation for such an array would need 
to maintain lists of voltage gain corrections for each an- 
tenna. This is straightforward and simply adds to the pro- 
cedure some bookkeeping as well as computational over- 
head. 

If large numbers of CLEAN components are involved, 
this technique is expensive in computer cycles and mem- 
ory. A simple expedient procedure is to "compress" the 
list by summing the flux density of all components derived 
from a given image grid cell and using at most one compo- 
nent per cell. An additional benefit of summing the com- 
ponents in each cell is the reduction of numerical round-off 
error in calculating the model response. 

3.2. Beam squint gain corrections 

This beam squint correcting algorithm has been imple- 
mented in the Obit package as task "Squint." The point- 
ing offset in antenna coordinates due to the beam squint is 
orthogonal to the line from the feed to the antenna center 
wi th its value gi ven by equation 2 using the values given 
bv lNapieJ (|l994h 

Thus, the offset for each polarization at the (E)VLA 
is given by 

, ,arcsecond. ,„, 

squint = ^237.56 x A( ) (3) 

TTieter 

where A is the wavelength (m). The value of squint is 
negative for LCP and positive for RCP. The offsets in RA 
and Dec are then: 



dx = squint * (— sm(feedAngle — x)) 
dy — squint * (cos (feed Angle — x)) 



(4) 



where feedAngle is the orientation of the feed on the feed 
circle, x is the parallactic angle and the sign of squint is 
appropriate for the polarization (jCotton fc Usonll2007a[ ). 
For a given CLEAN component, the distance from the 
pointing position corrected for beam squint is given by: 

d{x,y) ^ cos^^{sin{Soff)sin{dpnt) + 

cos{Soff)cos{Spnt)cos{aoff - apnt)} 

(5) 
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where apnt and Spnt are the RA and Dec of the antenna 
pointing, aoff = ace + x, Soff = Sec + V, and ace 
and See are the RA and Dec of the component. 

For frequencies above 1 GHz, the antenna pattern 
is approximated above the 5% level by a Jinc func- 



tiorlj which is in turn app roximated using §9.4.4 in 



Abramowitz fc StegunI ( 196 



P{x,y) = 4.0(0.5- 



n=l 



(6) 



where 



{1.496 X 10" 



25.0 



X d{x,y) X i^y 



and ci=-0. 56249985, C2=0.21093573, C3=-0. 03954289, 
C4=0.00443319, C5=-0.00031761, C6=0.00001109 and dant 
is the effective antenna diameter in meters with dgnt ^ 



24.5 m at a frequency of 1.4 GHz (jCondon et al.l 119981 ). 
The voltage gain for each antenna and polarization is then: 



voltage_gain ~ y^P{x,y) 



(7) 



As described above, the correction assumes a single 
azimuthally-symmetric antenna pattern describable by a 
1-D function. A single 2-D antenna pattern could be 
readily incorporated as a lookup table. Antenna-specific 
beam patterns require maintaining separate lists of com- 
ponent gains for each antenna. In the case of a large num- 
ber of components, this method can be time consuming. 
Therefore, in the implementation described here, a thresh- 
old is defined that specifies the cutoff between the levels 
that are processed using the high accuracy squint correc- 
tion model and those that are not. If the purpose of the 
application of this technique is solely to improve the dy- 
namic range in Stokes I by reducing the artifacts arising 
from the brightest emission in the field, it is not neces- 
sary to compute a highly accurate model for components 
corresponding to regions of weaker emission. By restrict- 
ing the accurate model calculation to the portion of the 
model representing the brightest emission and allowing 
faster, but less accurate calculations for the weaker emis- 
sion, the expense of this algorithm can be drastically re- 
duced. Because extended, low-level emission can produce a 
coherent error pattern it is important to keep this thresh- 
old rather low. In addition, if Stokes V imaging is desired, 
the beam squint correction must be applied to the entire 
Stokes I model. During CLEANing, whenever the peak 
in the initial residual image at the start of a major cy- 
cle exceeds the threshold value, the high accuracy model 
is required for all components. When the peak residual 
drops below this threshold, a lower accuracy model cal- 
culation is allowed if it is deemed to be faster. For self- 
calibration or other cases where the total CLEAN model 



^ The arbitrary cutoff at the 5% response is due to our lack 
of precise knowledge of the (E)VLA primary beam response 
at low levels; we are investigating ways to extend the squint 
correction to lower levels. 



is known, the components corresponding to a given image 
cell are combined and those whose flux densities exceed 
the threshold are processed with the high accuracy model, 
whereas those below the threshold may be processed with 
the faster method. In all cases, components corresponding 
to a given cell are combined before the model calculation. 
A nonzero threshold for using the more accurate calcula- 
tion might limit the dynamic range because of its effect 
on the amplitude self-calibration; a lower threshold leads 
to higher accuracy, although at the expense of execution 
time. 

3.3. Correction of individual antenna patterns 

Detailed differences among antennas will cause differences 
in their radiation patterns. For extremely high dynamic 
range, wide-field observations corrections for these dif- 
ferences among antennas will be necessary. This case is 
similar to the beam squint problem in that time-variable 
and directional gain corrections are needed. A straightfor- 
ward modification of the technique described in this pa- 
per would allow correcting for detailed antenna patterns, 
so that errors due to cross-polarization leakage could be 
corrected in the manner described here. 

3.4. Correction of antenna pointing errors 

As mentioned in the introduction, in the absence of point- 
ing errors, the effects of beam squint on RR and LL corre- 
lations will be approximately equal but of opposite sign. 
Thus, in this case, imaging of Stokes I (0.5 *(RR-|-LL)) will 
be largely, but not completely, insensitive to beam squint. 
However, in the presence of pointing errors, the magnitude 
of the beam squint on RR and LL will be different and will 
not be eliminated by imaging Stokes I, or by the technique 
described here without additional pointing corrections. 

At this writing, global, antenna-based pointing errors 
(i.e., coUimation errors) can be readily incorporated as our 
algorithm treats the beam-squint as a coupled (R,L) point- 
ing error. If post-observation corrections of time-variable 
antenna pointing errors could be determined, a straight- 
forward modification of this technique could be applied. 
Effective component gains for each antenna and polariza- 
tion could be kept, including the antenna pointing errors 
as well as the beam squint. A time and antenna depen- 
dent addition to the beam squint position offsets would 
be used in computing the effective component gains. This 
would simply require some bookkeeping and involve extra, 
straightforward computations. 

4. Verification using actual data 

4.1. Example 1: Off-axis observations of 3C84 

The beam squint algorithm was tested on VLA B configu- 
ration observations of 3C84 in which the source was alter- 
natively located on-axis, at the 1/2 power of the antenna 
beam (westward) and at the 1/2 power and 1/3 power of 
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Fig. 2. Time averaged 1.4 GHz Stokes V for observations witli 3C84 located at the half-power point westward of 
the pointing center. The left panel shows uncorrected data, the minimum is located where the Stokes V changes sign; 
positive and negative regions of Stokes V are indicated. The right panel shows the corresponding data after subtraction 
of the beam-squint corrected Stokes I model. The amplitude bias of the individual measurements is ~ 200 mJy. 



the antenna beam (southward). Observations were made 
in both the 1.4 and 5 GHz bands in spectral line mode, 
recording 15 x 390 kHz channels. The observations were 
made in July and August 2006. The data were amplitude, 
phase and bandpass calibrated in the normal fashion us- 
ing the on-axis pointings with 2-point interpolation of the 
gain solutions (amplitude, phase) to the other pointings 

The data were processed twice, once without and 
once with beam squint correction applied but otherwise 
processing was the sa me. All im aging used the "auto- 
center" technique (Cot ton fc Uso n 2007b), to improve 
the dynamic range. A single iteration of phase-only self- 
calibration with a solution interval of 30 seconds was 
followed by a single iteration of amplitude and phase 
self-calibration with a 2-minute solution interval. Self- 
calibration solutions were determined for an average of 
right- and left-hand circular polarization as without beam 
squint corrections, the independent right- and left-hand 
amplitude self-calibration would "correct" the data for 
3C84 at the expense of the rest of the field. In the fol- 
lowing paragraphs, we discuss the observations of 3C84 
at a frequency of 1.4 GHz with the source observed at 
the (westerly) half-power location of the antenna primary 
beam. 

A convenient way to see the effects of beam squint and 
to test the efficacy of the correction is to examine a plot of 
the time-averaged Stokes V amplitudes averaged as well 
over baselines and spectral channels (Figure [51 left panel). 
The instrumental polarization is fixed to the antenna and 
will rotate on the sky with parallactic angle causing a 
time- variable Stokes V signal in the averaged data. Since 
the source is not located at the phase center, this averaging 



must be done over amplitudes and therefore will suffer 
from a significant amplitude bias. Also, since the figure 
shows only the amplitude of the averaged visibilities, both 
positive and negative Stokes V appear with positive sign 
with a dip in the curve where the instrumental response 
changes the sign of Stokes V. 

Each data set was imaged using these corrections 
to derive a Stokes I model which was then subtracted 
from the data after performing the squint correction. A 
time-averaged Stokes V plot (Figure [2l right panel) will 
show the extent to which the procedure has worked. The 
200 mJy level in Figure [2] corresponds to the amplitude 
bias of the individual time and spectral samples. This fig- 
ure shows that the beam squint has been corrected to 
the level expected from the amplitude bias and within the 
scatter expected from the noise in the individual measure- 
ments. Similar results were obtained for the 5 GHz data. 

The Stokes I images are shown in Figure [3] with a 
square-root stretch of the lower pixel values in the image 
(peak 12.2 Jy/beam). The image made without beam- 
squint correction shows more residual artifacts than the 
image with the correction applied. To compare the two 
images quantitatively, the RMS noise was computed in 
identical boxes which are shown in Figure [31 The RMS in 
the image without correction was 1.47 mJy/beam whereas 
it was 1.16 m Jy/beam in the image made with the squint 
correction, a reduction of the level of artifacts by 21%. 

Figure m compares the Stokes V results. Since the data 
included observations over a range of parallactic angles, 
much of the instantaneous instrumental circular polariza- 
tion is washed out. Nevertheless, the image made without 
squint correction shows residual errors that are one order 
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Fig. 3. Left: Example 1. Stokes I in reverse gray-scale of the region around 3C84 imaged with no beam squint 
correction. The image is displayed with a square root transfer function and the scale is given by the wedge (in 
mJy/beam) at the top. The box gives the region over which the RMS was determined. 
Right: As left but with beam squint corrections. 





Fig. 4. Left: Example 1. Stokes V in reverse gray-scale of the region around 3C84 imaged with no beam squint 
correction shown with a linear display of -100 to 100 mJy/beam. The scale is given by the wedge at the top. 
Right: As left but with beam squint corrections. Notice the finer scale displaying the range of values -10 to 
10 mJy/beam. 



of magnitude larger than those in the squint-corrected im- 
age which is displayed using only the inner 10% of the 
pixel range of the uncorrected image. In the uncorrected 
image, the maximum and minimum values are well off the 
source but the maximum value in the corrected image is 



on-source. However, this analysis becomes confused as the 
core in this source is known to show circular polarization 
above the level shown in Figure [H Since the data were cal- 
ibrated using 3C84 and assuming no circular polarization, 
the polarization calibration is incorrect. Unfortunately 
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no other calibrators were included in the observations. 
The off-source noise levels computed using the same re- 
gion as in the Stokes I images are 7.9 mJy/beam and 
0.94 mJy/beam for the un-corrected and corrected images 
(respectively) . 

In this test, the beam squint correction leads to a sig- 
nificant improvement to the quality of the Stokes I image 
and to a very large improvement in the Stokes V image. 



4.2. Example 2: Imaging of full-track observations 

Our second example is a deep, full-track (two passes) 
1.4 GHz Hi observation of the field of the superthin 
galaxy IC2233 which contains significant background con- 
tinuum emission totaling '^ 1.1 Jy (observed, i.e. before 
correction for the attenuation of the primary beam) which 
is dominated by two bright "4C" sources (4C-I-46.17 at 
J081430.4-f455639.4 (henceforth the "upper-left" source) 
and 4C+45.15 at J081242.0+453651.3 (hereafter the 
"lower- right" source). The sources are observed with flux 
densities of 0.87 Jy and 0.18 Jy, uncorrected for the 80% 
and 34% (respectively) attenuation of the primary beam. 
These sources limit the dynamic range significantly. The 
data were obtained with the VLA in C configuration on 
2000 May 28-29 during a period of solar maximum. A full 
analysis of the spectral observation s has been published 
elsewhere ( Matthews fc UsonI 12008). We discuss here the 
imaging of the continuum emission with and without the 
squint correction as well as the limiting factors and resid- 
ual errors in the squint-corrected images. 



The observations were made in two overlapping fre- 
quency bands, each with 63 x 24.4 kHz channels, which 
allowed for 6 line-free spectral channels from the low- 
frequency side of the first (lower frequency) band and 14 
line-free channels from the high-frequency side of the sec- 
ond (higher frequency) band, after discarding 6 channels 
on each end that were contaminated by ghost images from 
the two strong "4 C" sources in the field (Uson .2007!) . 
The remaining line- free channels (6-1-14) were processed 
using conventional methods to make the images shown 
on the left-hand side of figures [5] and [S] Phase-only self- 
calibration was performed on a time scale of 20 seconds 
with the corrections applied to the data, followed by am- 
plitude and phase self-calibration performed on 5-minute 
time intervals after pre-averaging the RCP and LCP data 
with these corrections applied to both polarizations. The 
Stokes I image has an RMS noise level of 150 /iJy/beam 
and quite significant residual sidelobe patterns surround 
the two "4 C" sources. The variable squint results in a 
modulation of the apparent Stokes V signal as a func- 
tion of time. Thus, the uncorrected Stokes V image shown 
in the left-hand panel of figure [5] shows strong sidelobes 
that cannot be deconvolved as they correspond to a vari- 
able "source" that does not satisfy the average deconvolu- 
tion equation expected from the synthesized beam that 
corresponds to the full uv-coverage. The RMS noise is 



190 /iJy/beam but the noise distribution is clearly non- 
Gaussian. 

We have applied our squint-correcting procedure to 
these data. The procedure included automatic selection 
of cell-size and imaging grid of " fly's-eye" fields and out - 
liers, as well as "auto-centering" ( Cotton fc Usonll2007br ) . 
to improve the dynamic range. The procedure incorpo- 
rates phase-only self-calibration (20-second intervals) fol- 
lowed by amplitude and phase self-calibration (5-minute 
intervals). It is worth noting that the self-calibration op- 
erations used about one half of the number of degrees 
of freedom used with the conventional analysis because 
both bands were combined into one 86-channel "band" 
before applying our algorithm. The RMS noise in the re- 
sulting images is 110 /iJy/beam for the Stokes I image and 
105 /iJy/beam for the Stokes V image. 

The results are shown in the right-hand panels of fig- 
ures [5] and [6l Although the Stokes V image is largely con- 
sistent with noise with peak values at ±4.7(t, a faint pat- 
tern of residual sidelobes is visible surrounding the posi- 
tion of the upper- left "4 C" source. Some residual sidelobes 
to the second-strongest source remain in the Stokes I im- 
age as well. This lower-right source is located at a point of 
steep gradient in the primary beam response and is thus 
significantly affected by coUimation and pointing errors. 
We believe that such errors are likely to be responsible 
for our failure to correct this image in full. Of course, 
any uncorrected systematic effects will lead to errors in 
the self-calibration solutions. Despite the presence of such 
residual systematic errors, the Squint algorithm yielded 
an increase in dynamic range by more than one order of 
magnitude in Stokes V and a decrease of the rms noise in 
the Stokes I image of ~ 27%. 

5. Conclusions 

Taking advantage of the EVLA's increased sensitivity at 
lower frequencies (15 GHz and below) will generally re- 
quire imaging and deconvolving the entire primary beam, 
or at least those portions containing sources. Beam squint 
can produce artifacts well above the noise level. Other 
interferometers with off-axis feeds, including focal plane 
array feeds, will share this problem 

This paper describes a technique for correcting ra- 
dio intcrfcrometric observations for the effects of beam 
squint. The technique is a modification of the visibility- 
based CLEAN in which the Stokes I image is made in the 
usual way, but when the image model is subtracted from 
the data, corrections are made such that the model sub- 
tracted includes the effects of beam squint. Thus, any er- 
rors incurred in the initial imaging are corrected through 
iterative imaging and accurate model subtraction. Self- 
calibration is incorporated into the procedure in two steps 
with phase-only self-calibration performed first, subse- 
quently followed by amplitude and phase self-calibration. 

We have performed several tests of an implementation 
of this technique in the Obit task Squint. In all cases, 
careful initial calibration is critical. Observations of the 
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Fig. 5. Left: Example 2. Stokes I in reverse gray-scale of the field of IC 2233 imaged with no beam squint correction. 
Notice the deconvolution errors due to the two "4 C" radio sources. The image is displayed with the scale given by 
the wedge (mJy/beam) at the top. The rms noise is ^ 150 /iJy/beam. 

Right: As left but with beam squint corrections. The noise is reduced to ^ 110 /xJy/beam. Some artifacts remain 
that are centered on the "4 C" source on the lower-right side. 





Fig. 6. Left: Example 2. Stokes V in gray-scale of the field of IC 2233 imaged with no beam squint correction. Notice 
the artifacts due to the variable Stokes V "signal" at the position of the two "4 C" sources induced by the Earth's 
rotation of the (squinted) primary beam response with respect to the sky. This error pattern determines the RMS 
noise value of 190 /iJy/beam. The scale is given by the wedge (mJy/beam) at the top. 
Right: As left but with beam squint corrections. The rms noise is ^ 105 /iJy/beam. 



strong source 3C84 were made with the source at various 
locations in the antenna beam pattern resulting in strong 
instrumental circular polarization. Both 1.4 and 5 GHz 
test data sets were examined and we found that this tech- 
nique corrects the instrumental circular polarization due 



to the beam squint to a level that is limited by other ef- 
fects such as pointing errors and insufficient knowledge of 
the primary beam. 

We have presented detailed examples of imaging of sev- 
eral observations. Overall, applying the beam squint cor- 
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rection lowered residual artifacts in the vicinity of strong 
sources by more than 20% and led to very substantial im- 
provements to the Stokes V images. Our second example 
demonstrates the behavior of the algorithm on an observa- 
tion that requires high dynamic range. Here, artifacts re- 
main that can be ascribed to coUimation and pointing er- 
rors. It would be straightforward to extend this technique 
to correct for individual antenna patterns and known (i.e., 
determined elsewhere) antenna pointing errors. 

We are investigating the primary beam and pointing 
parameters of the (E)VLA antennas at several frequency 
bands and hope to achieve a more accurate characteri- 
zation of the squint correction which should lead to im- 
proved correction of the effects discussed here as well as of 
the next order of cross-polarization (leakage) corrections. 
We hope to report on those results in due course. 
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